home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / randist / flat.c < prev    next >
Encoding:
C/C++ Source or Header  |  2000-05-05  |  1.4 KB  |  54 lines
  1. /* randist/flat.c
  2.  * 
  3.  * Copyright (C) 1996, 1997, 1998, 1999, 2000 James Theiler, Brian Gough
  4.  * 
  5.  * This program is free software; you can redistribute it and/or modify
  6.  * it under the terms of the GNU General Public License as published by
  7.  * the Free Software Foundation; either version 2 of the License, or (at
  8.  * your option) any later version.
  9.  * 
  10.  * This program is distributed in the hope that it will be useful, but
  11.  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13.  * General Public License for more details.
  14.  * 
  15.  * You should have received a copy of the GNU General Public License
  16.  * along with this program; if not, write to the Free Software
  17.  * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  18.  */
  19.  
  20. #include <config.h>
  21. #include <math.h>
  22. #include <gsl/gsl_rng.h>
  23. #include <gsl/gsl_randist.h>
  24.  
  25. /* This is the uniform distribution in the range [a, b)
  26.  
  27.    p(x) dx = 1/(b-a) dx   if  a <= x < b
  28.    .....   = 0            otherwise 
  29.  
  30.  */
  31.  
  32. double
  33. gsl_ran_flat (const gsl_rng * r, const double a, const double b)
  34. {
  35.   double u = gsl_rng_uniform (r);
  36.  
  37.   /* A uniform distribution over [a,b] */
  38.  
  39.   return a * (1 - u) + b * u;
  40. }
  41.  
  42. double
  43. gsl_ran_flat_pdf (double x, const double a, const double b)
  44. {
  45.   if (x < b && x >= a)
  46.     {
  47.       return 1 / (b - a);
  48.     }
  49.   else
  50.     {
  51.       return 0;
  52.     }
  53. }
  54.